How perfect can graphene be? 
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We have identified the cyclotron resonance response of purest graphene ever investigated, which 
can be found in nature on the surface of bulk graphite, in form of decoupled layers from the substrate 
material. Probing such flakes with Landau level spectroscopy in the THz range at very low magnetic 
fields, we demonstrate a superior electronic quality of these ultra low density layers (no ~ 3 x 
10 9 cm -2 ) expressed by the carrier mobility in excess of 10 7 cm 2 /(V.s) or scattering time of r ~ 20 ps. 
These parameters set new and surprisingly high limits for intrinsic properties of graphene and 
represent an important challenge for further developments of current graphene technologies. 



Fabrication of graphene structures has triggered vast 
research efforts focused on the properties of two- 
dimensional systems with massless Dirac fermions. Nev- 
ertheless, further progress in exploring this quantum elec- 
trodynamics system in solid-state laboratories seems to 
be limited by insufficient electronic quality of manmade 
structures and the crucial question arises whether exist- 
ing technologies have reached their limits or major ad- 
vances are in principle possible. The substrate, and more 
general any surrounding medium, has been recently iden- 
tified as a dominant source of extrinsic scattering mech- 
anisms, which effectively degrade the electronic quality 
of currently available graphene samples [H, 0]. Despite 
significant advances in technology, including the fabrica- 
tion of suspended specimens 0,13], the realistic limits of 
the scattering time and mobility in graphene, achievable 
after elimination of major extrinsic scattering sources, 
remain an open issue. Experiments call for higher qual- 
ity samples, which are almost certainly crucial for possi- 
ble verification of interesting predictions concerning ba- 
sic phenomena of quantum electrodynamics (e.g. Zittcr- 
bewegung) or observation of the effects of interactions 
between Dirac fermions, (resulting e.g. in the appear- 
ance of the fractional quantum Hall effect). Even more 
simple effects such as lifting the degeneracy of the spin 
and/or pseudo-spin degree of freedom in very high mag- 
netic fields indicate the strong influence of the sample 
quality on the information that can be deduced from ex- 
periments [HQ. 

Recently, well-defined graphene flakes have been dis- 
covered in form of sheets, decoupled from, but still lo- 
cated on the surface of bulk graphite which naturally 
serves as a well-matched substrate for graphene [9(- In 
this Letter, we report on response from these flakes in mi- 
crowave magneto-absorption experiments and show that 
their Dirac-like electronic states are quantized into Lan- 
dau levels (LLs) in magnetic fields down to 1 mT and 
at elevated temperatures up to 50 K. The deduced un- 
precedented quality of the studied electronic system sets 
surprisingly high limits for the intrinsic scattering time 
and the corresponding carrier mobility in graphene. 



The cyclotron resonance has been measured in a high- 
frequency EPR setup in double-pass transmission config- 
uration, using the magnetic-field-modulation technique. 
A flake of natural graphite was placed in the variable tem- 
perature insert of the superconducting solenoid and via 
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FIG. 1: Magneto-absorption spectra (detected with field mod- 
ulation technique) of natural graphite specimen taken in the 
temperature interval 7.5 — 50 K at fixed microwave energy 
huj — 1.171 meV. The response of low electron concentration 
graphene layers decoupled from bulk graphite is seen within 
the yellow highlighted area. High CR harmonics of K point 
electrons in bulk graphite Q , with basic CR frequency ui c cor- 
responding to the effective mass m = 0.054mo 8] are shown 
by red arrows. The origin of features denoted by blue arrows 
is discussed in the text. 
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FIG. 2: Magneto-absorption spectrum (after removing a weak 
linear background seen in the raw data of Fig. [TJ measured 
at T — 25 K and microwave frequency tuu — 1.171 meV (a) 
in comparison with the LL fan chart (b), where the observed 
CR transitions are shown by arrows. The occupation of indi- 
vidual LLs is given by the Fermi-Dirac distribution plotted in 
the part (c). For simplicity, we considered only n-type dop- 
ing with Ef = 6.5 meV. The dashed lines show positions of 
resonances assuming c = 1.00 x 10 6 m.s -1 . 



quasi-optics exposed to the linearly polarized microwave 
radiation emitted by a Gun diode at frequencies 283.2 or 
345 GHz (1.171 or 1.427 meV). The absorbed radiation 
has been followed by either heterodyne detection (283.2 
GHz) or by the bolometer (345 GHz). To enhance the 
relatively weak response of graphene flakes, the modula- 
tion amplitude (AB ~ 0.5 mT) had to be chosen close 
to the CR width, which broadens resonances observed at 
lower magnetic fields. All spectra have been corrected 
for the remanent field of the magnet. 

Our main experimental finding is illustrated in 
Fig. [TJ The traces in this figure represent the magneto- 
absorption response of the natural graphite specimen at 
different temperatures, measured as a function of the 
magnetic field at fixed microwave frequency. They corre- 
spond to the derivative of the absorption strength with 
respect to the magnetic field since the field modulation 
technique has been applied. Strong at 7.5 K, but rapidly 
vanishing with temperature lines marked with red arrows 
in Fig. [T] can be easily recognized as cyclotron resonance 
(CR) harmonics of K point electrons in bulk graphite. A 
possible origin of the temperature dependent transition 
marked with the blue arrow will be discussed later on. 
The features of primary interest, which we argue are due 
to decoupled graphene sheets on the graphite surface, are 
seen at very low fields, within the yellow highlighted area 



of Fig. m 

The interpretation of these low-field data is schemat- 
ically illustrated in Fig. [5] The observed spectral lines 
(Fig. [5^) are assigned to cyclotron resonance transitions 
between adjacent LLs (|An| = 1) with energies: E n = 
s\gw(n)c^/ 2ehB\n\ — sign(n)£'i \J B\n\ 0, [ll[, charac- 
teristic of massless Dirac fermions in graphene sheets 
with an effective Fermi velocity c. This velocity is the 
only adjustable parameter required to match the ener- 
gies of the observed and calculated CR transitions. A 
best match is found for c = (1.00 ± 0.02) x 10 6 m.s -1 in 
fair agreement with values found in multilayer epitaxial 
graphene II 2. 1X311 or exfoliated graphene on Si/Si02 sub- 
strate 14j, [l5j, [l6(. As can be seen from Figs. [2}},c the 
multi-mode character of the measured spectra is directly 
related to thermal distribution of carriers among different 
LLs. The intensity of a given transition is proportional 
to the difference in thermal occupation of the involved 
LLs. Roughly speaking, the strongest transitions imply 
LLs in the vicinity of the Fermi level, which fixes at 
around 6-7 meV from the Dirac point. 

To reproduce the experimental data, we assume the 
absorption strength is proportional to the longitudinal 
conductivity of the system: 

Jn ~ fm 
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where /„ is occupation of the n-th LL and M m n — 
with a = 2 for 11 or m equal otherwise 
17J. The calculated traces in Fig. [3] have 
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been drawn taking 7 = 35 /LteV for the line broaden- 
ing, c = 1.00 x 10 6 m.s" 1 and E F = 6.5 meV. To di- 
rectly simulate the measured traces, the derivative of 
the absorption with respect to the magnetic field has 
been calculated taking account of the field modulation 
AB = 0.5 mT used in the experiment. In spite of its 
simplicity, our model is more than in a qualitative agree- 
ment with our experimental data, see Fig. [3] The cal- 
culation fairly well reproduces the experimental trends: 
Multi-mode character of the spectra, the intensity dis- 
tribution among the lines, as well as its evolution with 
temperature, and allows to estimate the characteristic 
broadening of the CR transitions. 

Our modeling could be further improved but at the ex- 
penses of additional complexity which we want to avoid 
here. Assuming magnetic-field and/or LL index depen- 
dence of the broadening parameter 7 and taking into ac- 
count the possible fluctuation of the Fermi level within 
the ensemble of probed flakes would improve the agree- 
ment between experiment and theory, particularly at low 
temperatures. Comparing both, the measured and sim- 
ulated traces, we are also more confident in the spec- 
tacular observation of the CR transition (involving the 
n = LL) at a magnetic field as low as 1 mT. Bear- 
ing in mind the small value of the extracted broaden- 
ing parameter one may conclude that LL quantization 
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FIG. 3: Background-removed magneto-absorption spectra 
from experiment at T — 25 and 40 K (black curves) in com- 
parison to the response at the microwave frequency hui = 
1.171 meV calculated for parameters c = 1.00 x 10 6 m.s" 1 , 
7 = 35 /xeV and Ef = 6.5 meV. The dashed lines show the 
calculated positions of the CR transitions. 



should survive in studied graphene layers down to the 
field of B = (j/Ei) 2 w 1 fiT. Hence, the magnetic field 
of the Earth of ~ 50 fiT is no longer negligibly small. 
Instead, it can open an energy gap at the Dirac point up 
to A k 0.3 meV, depending on the sample orientation. 

To crosscheck our interpretation, we have also mea- 
sured the spectra using a different (higher) microwave 
energy Hlu — 1.427 meV, see Fig. |4j Despite the weaker 
sensitivity of the experimental setup at this frequency, we 
can clearly identify the same set of inter-LL transitions 
simply shifted to higher magnetic fields. 

The finite Fermi energy Ep « 6—7 meV, corresponding 
to a carrier density of no rj 3 x 10 9 cm~ 2 , indicates that 
the probed layers are in thermodynamical contact with 
the surrounding material, which supplies these carriers. 
On the other hand, we find no signs of electrical coupling 
of these graphene layers to bulk graphite. Our experi- 
ments show that any possible energy gap opened due to 
this interaction at the Dirac point cannot exceed a few 
hundred fieV. The absence of this gap convincingly con- 
firms that we are indeed dealing with decoupled graphene 
and not with the H point of bulk graphite, where Dirac- 
like fermions are also present [lU but a (pseudo)gap of a 
few meV is expected 1|| and indeed observed |2(j. The 
temperature evolution of the measured spectra is another 
important indication allowing us to discriminate between 
the graphene and bulk graphite contributions. No (or 
very weak) temperature broadening of CR transitions is 
expected for graphene (3,[l3, 21 1, whereas the response of 
bulk graphite should follow the relatively strong decrease 
of the carrier scattering time, expressed by the average 
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FIG. 4: Positions of CRs at microwave energies huj — 1.171 
and 1.427 meV. Each position represents an average value ob- 
tained at different temperatures. One background-removed 
CR spectrum at each microwave frequency with a high num- 
ber of observed CRs is shown in the inset. The grey curves 
show positions of CRs for c = 1.00 x 10 6 m.s -1 . 



mobility, which reaches up to 10 6 cm 2 /(V.s) at low tem- 
peratures, but falls down by one order of magnitude at 
T « 50 K [§] . Indeed, this behaviour is observed in Fig. Q] 
Whereas the CR harmonics of K point electrons in bulk 
graphite seen in the spectra at \B\ > 20 mT, disap- 
pear very rapidly upon increasing T, the graphene-like 
features survive and their intensity is simply following 
the vanishing difference in the occupation between the 
adjacent LLs. It is worth noticing that the graphene- 
like signal, although always substantially weaker then 
the response from bulk graphite, has been observed for a 
number of different specimens of natural graphite. Me- 
chanical scratching of the sample surface and fast ther- 
mal cooling of the sample has been found to enhance the 
signal from decoupled graphene in comparison with bulk 
graphite, likely helping in detaching graphene sheets from 
the graphite crystal. 

Since the well-defined LL quantization in our graphene 
flakes is observed down to \B\ = 1 mT, see Fig. [2l we 
obtain via the semi-classical LL quantization condition 
[iB > 1 the carrier mobility /i > 10 7 cm 2 /(V.s), almost 
two orders of magnitude higher in comparison with sus- 
pended [1, 0, |i[ or epitaxial graphene [l3| . The LL broad- 
ening 7, obtained via comparison of our experiment with 
the simulated traces, allows us to estimate the scattering 
time t ~ 20 ps (t = ft/7), which significantly exceeds 
those reported in any kind of manmade graphene sam- 
ples, see e.g. Refs. 0,(1], and gives an independent esti- 
mation for the mobility /1 = ere 2 /Ep ~ 3x 10 7 cm 2 / (V.s) 
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in good agreement with the estimate above. Even though 
we cannot verify this estimate by a direct electrical mea- 
surement, a near correspondence of the scattering time 
derived from CR measurements and transport scatter- 
ing time was recently verified on samples with a sig- 
nificantly lower mobility 14} . Moreover, the estimated 



mobility should not decrease with temperature, as no 
broadening of CRs is observed up to T = 50 K, when 
CR intensities become comparable with the noise. This 
extremely high value of mobility combines two effects: 
the long scattering time r and a very small effective 
mass to = Ep/c 2 w 2 x 10~ 4 TOo-EF[meV]. Remarkably 
the same scattering time in a moderate density sample 
(no = 10 11 cm -2 ), would imply the mobility still remain- 
ing high, around fi«5x 10 6 cm 2 /(V.s), and compara- 
ble to best mobilities of two-dimensional electron gas in 
GaAs structures at these densities. 

The model we used here to describe the magneto- 
absorption response of Dirac fermions in graphene is 
amazingly simple, based on a one-particle approximation, 
and it is perhaps surprising that it is so well applicable to 
simulate the experimental data, particularly in context of 
the outstanding quality of the electronic system studied. 
At first sight, the observation of collective excitations, 
due to, e.g., magneto-plasmons could be expected in our 
experiments and only preliminary theoretical work ad- 
dresses the surprising approximate validity of Kohn's the- 
orem in graphene 22, 23[. Also size-confined magneto- 
plasmons are apparent in the microwave-absorption spec- 
tra of two-dimensional gas of massive electrons [13] j but 
apparently not seen in our experiments on graphene. 
This perhaps points out the qualitative difference in the 
plasmon behaviour in systems with quadratic and lin- 
ear dispersion relations 25|, 26[. On the other hand, we 
speculate that the transition, marked with blue arrows 
in Fig. Q] has a character of a collective excitation. It is 
characterized by a sensitivity to the magnetic field and by 
softening with temperature. The energy of this excitation 
can be deduced to scale roughly as B/T, and this points 
towards its magnetic origin. For example, B/T energy 
scaling is typical of a ferromagnetic resonance, see e.g. 



Ref. |27|. Nevertheless, the origin of the spectral feature 



marked by blue arrows in Fig. [T] remains an intriguing 
puzzle. We see this feature repeatedly on different sam- 
ples, but at slightly different energies, each time how- 
ever characteristically evolving with temperature. With 
its apparent magnetic characteristics one should seek its 
origin in impurities [28[, structural defects [2^] or edge 
states [Io[ in graphene, but possibly also on the surface 
of bulk graphite [31| . 

To conclude, graphene layers decoupled from bulk 
graphite have been probed in cyclotron resonance ex- 
periment, which offers an unambiguous evidence of ex- 
tremely high carrier mobility in graphene exceeding 
10 7 cm 2 /(V.s). This measurement significantly shifts the 
limits of intrinsic mobility in graphene [211 ] and poses 



a quest for further development in the technology of 
graphene fabrication. Graphene samples with mobilities 
comparable to the best GaAs samples [32| thus seem to 
be achievable. 
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